Supersymmetry and the positron excess in cosmic rays 



G.L. Kane"^, Lian-Tao Wang"^, James D. Wclls''^^ 
^'^^ Physics Department, University of Michigan, Ann Arbor, MI 48109 
''^'^ Phystcs Department, University of California, Davis CA 95616 
^'^'^ Lawrence Berkeley National Laboratory, Berkeley, CA 94720 



o 
o 

(N 

Oh' 
(D 



> 
oo 
cn 

00 

o 



Oh: 

I 

^ ■ 



X 



Recently the HEAT balloon experiment has confirmed an excess of high-energy positrons in cosmic 
rays. They could come from annihilation of dark matter in the galactic halo. We discuss expectations 
for the positron signal in cosmic rays from the lightest superpartner. The simplest interpretations 
are incompatible with the size and shape of the excess if the relic LSPs evolved from thermal 
equilibrium. Non-thermal histories can describe a sufficient positron rate. Reproducing the energy 
spectrum is more challenging, but perhaps possible. The resulting light superpartner spectrum is 
compatible with collider physics, the muon anomalous magnetic moment, Z-pole electroweak data, 
and other dark matter searches. 
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Recent experiments: Good solutions to the cosmolog- 
ical dark matter problem often involve hypothesizing a 
stable weakly interacting massive particle (WIMP). The 
particles populate galactic halos providing gravitational 
support to the unusual constant velocity profiles of many 
galaxies. 

Direct experiments continue to look for WIMPs scat- 
tering off nuclear target detectors. Indirect experiments 
rely on annihilation of ambient WIMPs that produce 
an excess above background of photons, anti-protons, 
positrons or neutrinos in cosmic rays. Each of these 
experiments has its unique experimental challenges, and 
its unique astrophysical assumptions and uncertainties. 
For example, to be successful the direct searches need a 
significant local density of WIMPs, whereas discovery of 
a monochromatic photon line from WIMP annihilations 
generally requires a cusping distribution near the galactic 
center. The charged particle signals {p and e"*") require an 
accurate model describing their propagation and energy 
loss from their source at WIMP annihilations in the 
galactic halo to the detector on earth. 

Our imperfect understanding of the dark matter dis- 
tribution and other astrophysics uncertainties makes it 
impossible to predict which signal would be the first 
to demonstrate evidence for WIMP dark matter. For 
this reason, all the different experiments designed for 
this purpose are interesting and necessary parts of a 
comprehensive search strategy. Once WIMPs are found 
all the experiments provide information about their 
properties and help to determine the WIMP relic density. 

Recently, the HEAT collaboration has found tantaliz- 
ing evidence for unexpected structure in the e+/(e^+e+) 
energy spectrum The first set of data from 

the 1994-1995 flights indicated a rise or bump in the 
positron fraction at energies above about 7GeV. Using 
a different instrument, with different systematics, the 
HEAT collaboration found in the data of their 2000 flight 



a similar rise. The consistency between the data sets adds 
further confidence in the measured energy distribution of 
the positron fraction. 

Attempts at a standard supersymmetry interpretation: 
One of the most compelling theories for WIMP dark 
matter is supersymmetry. R-parity conserving supersym- 
metry naturally provides a dark matter candidate in the 
lightest supersymmetric partner (LSP). 

In some models, such as "minimal supergravity" , the 
LSP is mostly bino (fermion superpartner to the hyper- 
charge gauge boson). For relatively light superpartners 
(mass near the weak scale) , one finds in large fractions of 
the parameter space of these models that a simple ther- 
mal history calculation will give an answer remarkably 
close to the Vlh? ~ 0.1 needed for an acceptable cold 
dark matter candidate. LSP annihilations into positrons 
can then be searched for in cosmic rays [^^. However, 
the standard supersymmetry model does not explain the 
HEAT data, for two reasons. 

First, the positron excess is most simply produced by 
LSP annihilations into W bosons ^-|^, one of which 
subsequently decays into a positron. However, binos do 
not couple to PF's and so this final state is suppressed 
compared to other final states. There is still the 
option of producing positrons from cascade decays of 
the other final states of bino annihilation. For example, 
annihilations into tau leptons can produce positrons from 
leptonic decays of t+ or from fragmentation of r jets. 
However, the total annihilation rate for binos is small. 
Although this is correlated with a reasonable Vlh^ , the 
annihilation rate is insufficient to produce a large flux of 
positrons to overcome expected backgrounds. Therefore, 
the positron fraction signal is not expected to be visible, 
unless we have underestimated important astrophysical 
parameters considerably. The HEAT data is likely not 
explained by bino LSP theories where the relic abundance 
of LSPs is accurately computed from a simple thermal 
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history of the universe. 

Higgsino and wino dark matter: What is needed to ex- 
plain the HEAT signal is a large relic abundance, a large 
annihilation rate, and a rising distribution of positron 
fraction at energies above about 7 GeV. Higgsino and 
wino LSPs may do this. They couple at full strength 
to the W boson and have a large annihilation rate. As 
long as they have mass above m\Y, higgsinos and winos 
will annihilate predominantly into WW final states and 
so can produce a large number of high-energy positrons 
from W + X. This has been discussed recently 

in the context of traditional supersymmetry models with 
large higgsino fraction LSP and anomaly mediation 
with wino LSP 0. 

Since the relic abundance correlates inversely with the 
strength of annihilation, there is still the worry that there 
will be too few of these LSPs in the galactic halo to an- 
nihilate with each other and produce a signal. However, 
that argument is based on a standard thermal history 
calculation which predicts I^lsp ^ 10"'^. Non-thermal 
sources and non-standard cosmologies have been found 
to produce a significant relic abundance independent of 
the thermal annihilation rate |p|-pT[ . It is one of the im- 
portant conclusions of this paper that the higgsino dark 
matter density probably must be understood outside the 
normal thermal evolution framework if the HEAT data 
is indicating LSP annihilations. 

Details of the positron signal: In order for the reader 
to understand our results we will briefly describe the 
assumed dark matter density profile we use to produce 
expected positron fiuxes from LSP annihilations. The 
dark matter halo is assumed to be spherically symmetric 
isothermal sphere whose density at a position r from the 
galactic center is 

p(r)=po4^ (1) 

where po = 0.3 GeV/cm'^ is the local LSP density, a = 3.5 
kpc is the core radius, and ro — 8.5 kpc is the distance 
of the earth from the galactic center. 

The flux i^g+ of positrons at the detector for higgsino 
or wino can be calculated from 

^ ^ ^/ deGME,e)Y^{av)jAiAe) (2) 
^ f 

where (crv) / is the annihilation rate of xx iiito the flnal 
state / = WW or / = ZZ, C?e+ {E, e) is the positron 
propagation Green's function, and (e) is the average 
positron energy distribution function for the final state / 
at the source (pre-propagation). The A^^ functions are 
normalized such that J de A^^ (e) is the average number 
of positrons in decays of the final state /. 
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FIG. 1. Positron fraction as a function of energy. The 
lower dashed line is the expected signal fraction with no LSP 
annihilations for a certain set of astrophysical assumptions 
described in the text. The solid lines also include the positrons 
and electrons from the annihilations in the galactic halo of 
LSPs with mass = 83 GeV and — 200 GeV, with boost 
factors of 2.7 and 3.9 respectively. The 1994-1995 HEAT data 
is represented by the solid line cross-hairs, and the 2000 HEAT 
data by the dashed line cross-hairs. 

We utihze DARKSUSY [|| for calculating the flux. 
We also tested the results by simulating the A^^ (e) from 
Pythia. The Green's function Gf,+ {E, e) can be extracted 
from refs. |^,[l^. The numerical values we used are 
from with energy loss time te = 10^^ sec, and with 

energy-dependent diffusion constant 

K{e) = 6.1 X 10^^ (yc^)°^ 

As indicated earlier, the thermal relic abundance is 
much too small to be of cosmological signiflcance, but 
non-standard mechanisms can save the higgsino and 
wino as dark matter candidates. From here on we 
assume that the local density po is made up entirely of 
neutralino dark matter, assuming a non-thermal source 
for the LSPs such as from late decays of very heavy 
gravitinos [^. Then we no longer need to concern 
ourselves with neutralino relic abundance, since the value 
of Po captures all the information we need about LSP 
abundance in our positron flux calculation. An obvious 
consequence of this approach is that we do not rescale po 
according to the thermal relic abundance calculation. 

Fig. shows the positron fraction energy distribution 
for the HEAT data |^,^, expected distribution with no 
LSP annihilations, and expected distribution with LSP 
annihilations subject to the above assumptions. The plot 
is made for higgsino/gaugino mixed scenario with LSP 
masses of — 83 GeV and 200 GeV. 

We have normalized the positron distribution to 0.06 
at 10 GeV for each value of the LSP mass. To do this 
we had to arbitrarily multiply the flux calculated with 
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FIG. 2. The solid line is the average positron energy distri- 
bution from cascade decays resulting from LSP annihilations 
into W-boson pairs. The mass of the LSP is 100 GeV in this 
example. The dashed line tracks the positrons from — > 
direct decays, and the dotted line from — > fi'^ /t'^ — > 
direct decays. 



the above assumptions by a factor ( "astrophysical boost 
factor") of 2.7 and 3.9 for = 83 GeV and 200 GeV, 
respectively. An extra boost factor less than about 10 
is probably well within the astrophysical uncertainties 
of parameters used to calculate the flux. This gives us 
confidence that low-mass LSPs mw < ra^ < 200 GeV are 
worth pursuing as possible interpretations of the positron 
energy distribution and the total flux. 

The HEAT data appears to show a dip in the positron 
energy fraction near = 7 GeV. A dip would indicate 
that a signal should have a large bump in its positron 
distribution at energy above 7 GeV. XX ~^ W^W" 
annihilations are the best hope to produce a bump in the 
positron spectrum from LSP annihilations, since 
e^v decays lead to a peak in the positron spectrum at 
high energies. However, there are numerous other sources 
for positrons in W decays, including cascades from r 
and /i leptons, and decays of pions in jets. In Fig. || 
we show the average positron energy distribution from 
XX ~^ W'^W~ annihilations, simulated using Pythia 
results from e+e^ W'^W~ . The lack of a peak 
in this distribution clearly indicates that simple LSP 
annihilations cannot reproduce a strong peak in the 
positron energy spectrum. 

Therefore, in order for the LSP annihilations to be 
consistent with the data we have to assume that there 
is no significant dip in the data, but rather a change in 
slope. Given the error bars for the HEAT data points, 
this possibility is not out of the question. In this case, the 
signal arises from an LSP-induced positron distribution 
that is somewhat flatter than the background positrons. 
The prediction is best fit to the data when the number 
of signal positrons starts to become a significant fraction 
of the total positron rate at > 5 GeV. 

We briefly mention here another interpretation of the 



data which is somewhat fine-tuned, but would be more 
consistent with strong peaking in the positron energy 
spectrum. The electron sneutrino is stable or nearly 
stable if its mass is extremely close to the LSP's. In this 
case one could imagine x^* ~^ W~e^ annihilations with 
the positron energy peaked at i?e+ — m^(l — m^/4TO^). 
Numerically, to get a sharp peak at about 8 GeV requires 
TOj> -I- « my/ + 10 GeV. Whether the neutralino or 
sneutrino is the lightest would not be important. We 
have not carefully studied this possibility, although we 
recognize that the LEPII collider data would severely 
constrain it, and maybe even rule it out. 

Correlated phenomena: In the previous sections we 
have concluded that 

• Traditional supersymmetry with thermal relic 
abundance near Q.h? ~ 0.1 cannot yield an excess 
of positrons above background because the annihi- 
lation rate is too low, 

• A higgsino or wino LSP with mass mw < 

< 200 GeV could yield a significant excess of 
positrons above background provided the relic 
abundance is from a non-thermal source, 

• A generic higgsino or wino interpretation is consis- 
tent with the HEAT data only if no strong dip is 
present, 

• Uncertainties in the data and in astrophysical 
processes such as positron production, propagation 
and modulation means we may not need a new- 
physics interpretation of the HEAT data, although 
our current understanding suggests we do. 

If the higgsino or wino interpretation of HEAT data is 
correct, we should expect other correlating phenomena 
that can be measured and quantified. First, it is well- 
known by now that the recent excess [l4| in the muon 
anomalous magnetic moment is consistent with light 
supersymmetry ]l5|-p^. Particularly, a light higgsino 
and large tan /3 are helpful to get a large supersymmetric 
correction, since the higgsino-smuon-muon vertex is a 
tan/3 enhanced chirality flip. The supersymmetric con- 
tribution to (5a^ with higgsino or wino LSP is 



Sa^ < 14 tan /3 



100 GeV y 



X 10" 



(4) 



Equality in eq. (|j) is attained when all other sparticle 
masses are very close to the LSP mass. The measurement 
minus the Standard Model contribution is Sa^, = (41 ± 
16) X 10"^". Therefore, moderate tan/3 and low mass 
superpartners have no difficulty recovering the central 
value for the measured Ja^. 

Colliders can also search for higgsino and wino cold 
dark matter. These searches are notoriously difficult 
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because there is no guarantee that visible superpartners 
have mass close to the LSP mass and are therefore 
accessible by the colliders. In the case of higgsinos and 
winos, there are charged particles nearby, or W~^; 
however, they are almost degenerate in mass to the LSP. 
The production of e+e^ H^H^ may be high, but the 
final state of two soft pions from H^n^ is very 

difhcult to find. Searches have been conducted, and the 
mass limits for these sparticles are about mw just 
below the interesting region for the HEAT signal. Future 
lepton colliders will have a much higher mass reach, and 
hadron colliders will be useful if other superpartners are 
produced ||l9|] . 

Other astrophysics experiments may also see evidence 
for wino or higgsino cold dark matter. Nuclear target 
detectors have limits that are already sensitive to wino 
or higgsino LSPs in some parts of parameter space . 
However, these limits depend on squark masses, heavy 
Higgs boson masses, etc. which feed into the spin- 
independent nucleon-LSP scattering cross-section, and 
which have little to do with the positron fraction pre- 
diction. Therefore, it is difficult to predict how sensitive 
next generation cryogenic detectors will be to light 
higgsinos and winos. A possible signal not in conflict 
with our results has been reported at DAMA pl] , |2^ . 

On the other hand, loop-induced annihilations of 
XX ~* 77 very high for higgsinos and winos. One 
therefore expects a monochromatic photon signal to 
arise from annihilations of winos and higgsinos in the 
galactic halo As mentioned at the beginning, 

the astrophysical uncertainties of this calculation are 
quite different than the positron fraction calculation. 
Therefore, it is difficult to predict if experiments such 
as GLAST will see a signal, but we do expect so if the 
HEAT results are due to LSP annihilations. 

Finally, it has been suggested recently that the pre- 
cision electroweak data is more consistent with light 
superpartners p5| . Sncutrino masses below mw and 
slepton masses just above the experimentally allowed 
region are the most important requirements for the 
successful fit to data. Light gauginos are also helpful, but 
not as critical. Therefore our interpretation of the HEAT 
data, which requires superpartner masses near miy, is 
not only compatible with the precision electroweak data, 
but may be encouraged by it. This is another reason 
why the LSP interpretation is worthwhile pursuing even 
though it has difhculty reproducing the precise structure 
of the data. 
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